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Abstract 

We studied the transition between the antiferromagnetic and the surface spin-flop phases of a 
uniaxial antiferromagnetic [Fe(14 A)/Cr(ll A] x 2o superlattice. For external fields applied parallel 
to the in-plane easy axis, the layer-by-layer configuration, calculated in the framework of a mean- 
field one-dimensional model, was benchmarked against published polarized neutron reflectivity 
data. For an in-plane field H applied at an angle tjj ^ with the easy axis, magnetometry shows 
that the magnetization M vanishes at H = 0, then increases slowly with increasing H. At a 
critical value of H, a finite jump in M(H) is observed for if) < 5°, while a smooth increase of M 
vs H is found for if) > 5°. A dramatic increase in the full width at half maximum of the magnetic 
susceptibility is observed for ip > 5° . The phase diagram obtained from micromagnetic calculations 
displays a first-order transition to a surface spin-flop phase for low ip values, while the transition 
becomes continuous for ip greater than a critical angle, VWx ~ 4.75°. This is in fair agreement 
with the experimentally observed results. 

PACS numbers: 75.70.-i, 75.50.Ee, 75.10.-b, 75.30.Kz 
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I. INTRODUCTION 



It is well known that when a magnetic field applied along the easy axis of a uniaxial 
antiferromagnet exceeds a critical value H BSF = ^2H e Ha + H\, where H E is the exchange 
field and Ha is the anisotropy field, the system undergoes a first-order phase transition to a 
bulk spin-flop phase, characterized by sublattice magnetizations nearly perpendicular to the 
field direction.-^ In the case of a uniaxial antiferromagnet with one or two surfaces, which 
break the translational invariance in the direction perpendicular to the surface plane, the 
problem of determining the ground-state spin configuration in the presence of an external 
magnetic field applied along the in-plane easy axis was theoretically posed a few decades 
ago. The first model to be investigated was that of a semi-infinite stack of ferromagnetic 
planes, antiferromagnetically coupled and subject to a magnetic field antiparallel to the 
magnetization of the surface plane. For this system, when the ratio r between Ha and H E 
is very small (r = Ha/H e -C 1), the onset of a surface spin-flop phase was predicted using a 
continuum approximation.- 1 ^ This phase is characterized by a canting near the surface and 
is stable for field values H greater than a critical value Hssf = J H e Ha + H\ fa Hbsf/ V%- 
Nearly a decade ago, the existence of such a surface spin-flop phase was criticized^ because 
a discrete, nonlinear map approach showed that the instability of the antiferromagnetic 
configuration at Hssf simply leads to an interchange of the two sublattices.- 1 ^ Subsequently 
Pokrovsky and Sinitsyn 7 showed that, for a semi- infinite film with r<l,a quite similar 
result can be obtained also in the continuum approximation, provided that appropriate 
boundary conditions are assumed. 

The case of a finite stack with an even number N of planes (but still r <C 1) is quite 
different since, for Hssf < H < Hbsf, the system tends to realize an inhomogeneous 
configuration with the magnetizations of both surfaces parallel to the field direction, and a 
domain wall is thus located in the center of the stack. For H > H E sf the stable state of the 
finite film is a different inhomogeneous configuration, with the inner spins assuming nearly 
the bulk spin-flop configuration and the surface ones less deviated with respect to the field 
direction.- 

Fullerton et a/. 8 showed that the Fe/Cr(211) superlattice, obtained by the repetition of N 
ferromagnetic iron layers antiferromagnetically coupled through the Cr spacer, constitutes an 
excellent model system of a finite uniaxial antiferromagnetic film. Since then, a great number 
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of paper^^^^^^^i&i^^^ have been devoted to the study of the surface spin- flop phase 
transition. In fact, for sufficiently low thickness of the Fe layers, the system possesses a 
dominant uniaxial in-plane anisotropy along the Fe[ -1 1] direction, with Ha of the same 
order of magnitude as He- For example, the Fe/ Cr(211) superlattice^ with thickness tp e = 40 
A and tci = 11 A had an anisotropy-to exchange ratio r = Ha/He w 1/4, while for the 



system investigated in the present work and in Ref. [15|, with tp e = 14 A and tci = H A, 
one has r « 1/10. This is a major difference with respect to ordinary antiferromagnets, 
like MnF 2 , where r is usually much smaller (r rs 1/100). The consequence of an increased 
value of r in superlattices with respect to bulk antiferromagnets was investigated both 
experimentally- and theoretically^^- in Fe/Cr(211) superlattices with r »s 1/4, and was 
found to introduce a sequence of sudden jumps in the field dependence of the magnetization, 
in addition to the surface spin-flop jump occurring at H$sf- 

Recently, an accurate and systematic study of the phase diagram of a uniaxial antifer- 
romagnetic film with an even number of planes was performed by Rofiler and Bogdanov 19 
using an efficient conjugate gradient minimization technique in the case of the external 
magnetic field applied precisely along the easy axis. For r < 1, they found that there 
is only a first-order transition from the collinear antiferromagnetic (AF) phase to a sym- 
metric, inhomogeneously flopped phase with the spin-flop (SF) located in the center of the 
film. For r < 1, their calculations of the spin configuration confirmed previous theoretical 
findings. In fact they found a series of canted, asymmetric phases (Q), separated by 
first-order transition lines, between the AF and the SF phases.~ Within these intermediate 
C phases, the ground state of the system evolves from a canted configuration with a flop 
localized near one of the surfaces (Ci) to other configurations (C 2 , C 3 , ...) where the flop 
moves into the center, causing abrupt variations of the magnetization as the field intensity 
is increased.— 1 ^ Finally, upon further increasing the anisotropy (r > 1), they found that 
only first-order transitions between collinear (antiferro-, ferri- and ferromagnetic) states are 
possible. 19 

While the magnetic phase diagram of the finite AF film has been extensively investigated 
both theoretically and experimentally in the case of H applied parallel to the easy axis, 
no such studies are known for the case for which H is applied in-plane along an arbitrary 
direction with the easy axis. In the bulk case, the field-induced phase transition of a uniaxial 
antiferromagnet in the presence of a skew field forming an angle ip ^ with the easy axis 
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was theoretically studied by Rohrer and Thomas 20 using a mean-field approach and then 
by Fisher et al.—^ Neglecting zero-point motion effects, they determined the equilibrium 
configurations as a function of the skew field, and found that the phase boundary between the 
AF and the bulk SF phase extends only to a maximum angle with respect to the easy axis, 
where it ends in critical points. More precisely, they predicted the first-order SF transition 
to become continuous (i.e., second-order) for ip greater than a critical angle ^ max (bulk) = 
tan -1 [H a/ '(^He — Ha)]- This expression, first developed for a tetragonal system, was later 
found to be valid for an orthorhombic system, when ip is restricted to a plane comprising 
the easy and intermediate axes.— 

The small values of ?/>max(bulk) in ordinary bulk antiferro magnets (amounting to a few 
tenths of a degree) made the observation of the crossover from first- to second-order in the 
transition difficult. Early results on M11F2 were, in this respect, only qualitative. More 
direct evidence of the crossover character of the transition was provided by Butera et alM 
by measuring the staggered magnetization of of MnCl2-4H 2 0. However, the existence and 
nature of a bicritical point was proven by measuring the critical magnetic scattering of 
a number of systems, notably CuCl2-2H20, 25 CsMnBr3-2H20,— and the above mentioned 
MnCl2-4H 2 0. On intuitive grounds, one might expect that in Fe/Cr(211) superlattices, a 
similar crossover effect in the order of the surface SF transition should be present and should 
be more easily observable owing the higher value of the ratio r = Ha/ He (r ~ 1/10 or more) 
compared to ordinary bulk antiferromagnets like MnF 2 , where r « 1/100. 

The aim of the present paper is twofold: i) to extend to the film (i.e. finite) case the 
theoretical study of the magnetic phase diagram of a uniaxial antiferromagnet in a skew field, 
and ii) to provide experimental evidence that, in contrast with bulk antiferromagnets, in 
the case of the [Fe(14 A/Cr(ll A)] X 2o superlattice, the crossover in the order of the surface 
SF transition with increasing skew field might be observed. Concerning the theoretical 
analysis, in Sec. II we present two different approaches, the nonlinear map method and the 
Landau-Lifshitz-Gilbert micromagnetic simulation, for the ground state and full hysteresis 
calculations, respectively. In Sec. Ill experimental results are presented and discussed, 
obtained by different techniques, both for H parallel to the easy axis (ip = 0) and for a skew 
field (ip ^ 0). Finally, conclusions are drawn in Sec. IV. 
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II. FILM MODEL AND THEORETICAL FRAMEWORK 



We consider a film made of an even number, N, of parallel ferromagnetic planes that are 
antiferromagnetically coupled one to each other by a nearest neighbor exchange interaction 
and subject to a uniaxial in-plane anisotropy. For thin magnetic layers, the magnetostatic 
dipolar interaction is known to confine the spins to the film plane, so that, at equilibrium, the 
spins are necessarily in-plane and the dipolar energy is zero. One can therefore characterize 
the spin configuration of the z-th plane by only one parameter, namely the angle fa that the 
magnetization of the z-th plane forms with the z axis (i.e., with the easy anisotropy axis) 
within the film plane xz. A magnetic field H is applied in the film plane along a direction 
that forms an angle ip with z. At T = K the energy density of the system takes the form: 

e = E/{g^ B SN\\) 

N 

= YjI H e cos(0j - 0j_a) - H A cos 2 fa 

i=l 

- 2H Z cos fa - 2H X sin fa] , (1) 

where N\\ denotes the number of spins within each film layer, and Hz = Hcosip, and 
Hx = H sin ip. 

The equilibrium spin configurations can be obtained from ([1]) by ^-derivation (i = 
AO 




= - H E sm(4>i - 0i_i)(l - Sij) 



- H E sin((p i+1 - (pi){l - S N)i ) 
+ 2Ha sin (pi cos (pi 

+ 2H Z sin fa - 2H X cos fa (2) 

In the present work, the ground state of the one-dimensional (ID) model ([1]) that approx- 
imates the film has been determined using two different theoretical methods, namely: i) an 
integration of the Landau-Lifshitz equation for the spin chain, introducing a damping coef- 
ficient in order to reach the magnetic ground state (see Sec. II. A), and ii) a reformulation 
of Eq. (J2j), which provide the equilibrium conditions in terms of a nonlinear map with op- 
portune boundary conditions at the film surfaces (see Sec. II. B). A comparative discussion 
of the two methods is made in Sec. II. C. 
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A. Spin configuration via integration of the Landau-Lifshitz-Gilbert equation 

The Landau-Lifshitz-Gilbert (LLG) equation 27 



describes the physical path of the magnetic moment M in a field H e g- Here 7 and «g 
are the gyromagnetic ratio of the free electron spin and the Gilbert damping coefficient 
respectively, M s is the saturation magnetization, and H e fj = —(1/M s )(de/&M.) is the local 
effective magnetic field. For thin magnetic layers, a large value of ac is appropriate be- 
cause the demagnetizing field confines the magnetic moment to the film plane, suppressing 
the gyromagnetic precession and leaving an in-plane rotation of the moment towards the 
direction of the local effective field. Thus, the integration of the LLG equation becomes 
a simple relaxation of the magnetic moment along the energy gradient. In the numerical 
calculations, we iteratively rotate spins that represent individual Fe layers by an amount 
proportional to the torque M x (M x H c ff) at each instant. H e ff is then evaluated from the 
resulting configuration and applied to the next iteration. Upon reaching convergence, the 
stability of the equilibrium is tested by evaluating the eigenvalues of the stability matrix 
Mij = d 2 e/(d<frid<f)j)M^ All eigenvalues of the stability matrix must be positive for the state 
to be stable. In the event of an unstable equilibrium, the configuration {0,} is displaced by 
a random, small fraction along the eigenvector direction for which the eigenvalue is nega- 
tive, and the relaxation process starts anew until the system reaches a stable local energy 
minimum. 

B. Determination of the energy minima via the nonlinear map method 

A different approach for the determination of the ground state of the magnetic system 
described by Eq. (Q]) was proposed some years ago.-^ 1 ^ It is based upon a reformulation of 
the equilibrium conditions (T2J) of the magnetic film model in terms of a discrete nonlinear 
map, where the film surfaces are introduced via opportune boundary conditions.— 

We start by introducing the new variable Sj = sin(0j — <f>i-i), so that the conditions for an 
equilibrium spin configuration, Eq. (j2J), can be rewritten as a 2D nonlinear recursive map. 




(3) 
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For 1 < i < N one has 



s i+1 = Si - (H A /H E ) sin(20i) 
- 2(H z /H E )sm& 
+ 2(H X / H E ) cosfa, 

<f>i+i = <j>% + sin~ 1 (sj + i). (4) 

The trajectories in (<ft, s) space are associated with equilibrium configurations, while the 
fixed points of the map correspond to uniform ground states of the infinite system and are 
second-order [i.e., they satisfy the relation (<p n+2 , s n+2 ) = (4>ni s n)] owing to its AF nature. 
We denote by (at, 0) the ground state configuration of the infinite system in the presence 
of a field of arbitrary direction. The angles a and (3 that the magnetizations of the two 
sublattices form with the easy axis can be determined by numerically solving the following 
problem of extremum in 2D: 

dej da = = H E sin(/3 - a) + (H A /2) sin(2a) 
+ Hz sin a — Hx cos a, 

de/d(3= = H E sin(a - (3) + (H A /2) sin(2/3) 

+ H z sin (3 - H x cos /3. (5) 

From Eq. ([5]) it is readily found that the second-order fixed points of the map are 
FPi=(/5, sin(/5 — a)) and FP2=(a, — sin(/5 — a)). In order to study the map behavior in 
the proximity of the fixed points, it is useful to perform a linear stability analysis of the 
doubly iterated map. It is worth noticing that energetically stable configurations (i.e., with 
a positive definite Hessian) are associated with topologically unstable (i.e., hyperbolic) tra- 
jectories in phase space. And, vice versa, energetically unstable configurations are associated 
with topologically stable (i.e., elliptic) trajectories in phase space.— 

At this point, the presence of the two surface planes, signaled by the terms with the 
Kronecker's 5's in Eq. (T5]), is taken into account via opportune boundary conditions: i.e., 
we introduce two fictitious planes i = and i = N + 1, so that we can assume the bulk 
equations, Eq. (Tj0), to be valid even for the surface planes i = 1 and i = N, provided that 
the following equations are satisfied: 

si = sin(0i - (f) ) = 
7 



s N = sm((f) N+1 - <f) N ) = 0. (6) 

In this way, among all trajectories obtained from the map equations (jlj), only those satis- 
fying the boundary conditions ([6]) represent equilibrium configurations of the film with a 
finite number iV of planes. In practice, the physical trajectories of the film must have two 
intersections with the s = line, separated by exactly N steps of the recursive mapping.-&ii 
Using this nonlinear map method, one is able to numerically determine all the stationary 
configurations of the film very rapidly and within machine precision. To find the ground 
state among the various calculated equilibrium configurations, it is necessary first to perform 
a linear stability analysis of the obtained solutions through the evaluation of the Hessian 
and then to compare the energies of the different metastable states in order to choose the 
lowest energy one. 

C. Comparison of the two methods 

By definition, the two theoretical methods described in Sec. II. A and II. B give the same 
results for the field evolution of the ground state from the AF to the SF phase. For the 
sake of generality, a comparison between the two methods and a critical discussion of their 
respective advantages and drawbacks will now be presented. 

The nonlinear map method allows one to determine all the equilibrium configurations 
of the film at zero temperature. Thus, comparing the energies of the various equilibrium 
configurations, one can easily determine the ground state. However, in the presence of two 
or more energetically equivalent ground states of different configurations, their knowledge 
is not enough to predict when a given equilibrium state is abandoned in favor of another 
equilibrium state as the intensity of the external field is varied. Such a process cannot be 
simulated in the framework of the nonlinear map method, whereas, in the approach based 
on the integration of the Landau-Lifshitz-Gilbert equation, the inclusion of the damping 
coefficient allows the system to evolve towards a stable local energy minimum. It should be 
noted that, in the framework of the latter method, the magnitude of perturbation applied 
during the stability test (i.e. the fractional displacement of the configuration {(pi} in the 
event of an unstable equilibrium) may be crucial in determining which one of these states is 
eventually reached. 
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To this regard, we note that the problem of the determination of the actual state of the 
system may be nontrivial in the case of uniaxial anisotropy field and exchange field with 
comparable intensities, r = Ha/H e 1, as in Fe/Cr(211) superlattices for a suitable choice 
of the layer thicknesses.- In fact the metastability region turns out to be strongly amplified 
with respect to the case r<l, and many metastable states with slightly different energies 
may be present, so that a tiny increment in the value of H may cause an abrupt change in 
the ground-state configuration.—^^ Such a peculiar feature is usually signaled by a chaotic 
aspect of the nonlinear map in the ((f), s) phase space.- 1 ^ In the extreme case of r ^> 1, 
fractal structures were predicted to appear both in the distribution of magnetic moments 
and in the energy spectrum.— 

As regards the [Fe(14 A)/Cr(ll A)] x2 o superlattice under study, in Section III it will 
be shown that this system is characterized by a moderate value of the anisotropy (r « 
1/10). As a consequence, the metastable state to which the system eventually relaxes in 
the micromagnetic calculations is not very sensitive to the magnitude of the perturbation 
applied during the stability test, so long as it is reasonably small. Thus, owing to the low 
value of r, the two theoretical methods were found to give quite similar results for the ground 
state of the Fe/Cr superlattice. (For the sake of precision, the results reported in Fig. 1 were 
obtained by the approach based upon the integration of the LLG equation, while those in 
Fig. 2, Fig. 6 and Fig. 7 were obtained via the nonlinear map method.) 

Finally, it is worth observing that the theoretical results refer to the ground state, while 
experiments (as described in Section III) are performed at room temperature T am b. The 
influence of a finite temperature on the spin-flop transition of a classical, simple cubic lat- 
tice, uniaxially anisotropic Heisenberg antiferromagnet was investigated some time ago.~ 
Roughly speaking, comparing the T = results with the experimental ones at T am b is rea- 
sonable as long as T am b ^5 T&, the bicritical point where the AF and SF ordered phases meet 
in the (H, T) phase diagram. For the square lattice, recent work 30 indicated that a very 
narrow disordered phase may intervene between the AF and the SF phase down to quite low 
temperatures, leading to the definition of a tetracritical point T t > XJ,. In both the D = 3 
and D = 2 cases, the multicritical point in the (H, T) phase diagram was found to be a sub- 
stantial fraction of T/v, the AF-paramagnetic transition temperature for H = 0. The Fe/Cr 
superlattice under study cannot be described by either of the two afore-mentioned models, 
since it is made of strongly ferromagnetic iron films, antiferromagnetically coupled through 
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the chromium spacer. However, denoting by T mc the multicritical point in the (H, T) phase 
diagram, one can expect T mc < T/v also in the superlattice case. In fact, on the basis of 
a mean field theory estimate, T/v is expected to be much greater than T amb for the Fe/Cr 
superlattice, so that one can guess also the condition T amb < T mc to be satisfied. 

III. EXPERIMENTAL RESULTS 
A. Sample characterization 

The preparation and characterization of epitaxial Fe/Cr superlattices are similar to those 
described in Ref. jsl The [Fe(14 A)/Cr(ll A)] x2 o superlattice was prepared^ by dc magnetron 
sputtering onto a single- crystal MgO(llO) substrate. To assure epitaxy with the substrate, 
a 200- A buffer layer of Cr was first deposited at 400 °C, then the superlattice was deposited 
at 100 °C and found to grow with a (211) orientation. Finally, a 100-A capping layer of Cr 
was deposited to protect the sample. The epitaxy and the smoothness of the superlattice 
were checked by x-ray diffraction and found to have an interfacial roughness of ~ 4 A. 
Extensive magnetic characterizations were performed by means of magnetometry, as well as 
by magnetoresistance measurements. For Fe film thickness t Fe = 14 A, a strong, in-plane 
surface anisotropy Kg = 0.06 ergs/cm 2 was found to develop along the [0 -1 1] direction, 
leading to a uniaxial in-plane anisotropy Kjj = 2Ks/t Fe = 8.6 10 5 ergs/cm 3 (compared 
to the bulk crystalline anisotropy K x = 4.7 10 5 ergs/cm 3 ). Using the value M s = 1740 
emu/cm 3 of bulk Fe at T = K, one obtains = 2Kjj/M s = 985.2 Oe. The coupling 
between ferromagnetic layers was found to oscillate as a function of the thickness of the Cr 
interlayer. For to = 11 A an AF exchange coupling with strength Jaf = —1.194 ergs/cm 2 
was estimated,^ leading to H E = 2\J AF \/(t Fe M s ) = 9802.9 Oe. 

For the [Fe(14 A)/Cr(ll A)] x2 o superlattice under study, previous SQUID measurements^ 
of M(H) vs H, applied in-plane along the easy axis of the sample, showed that M is zero 
for zero field. As H was increased, the instability of the AF phase was signaled by a steep 
increase of M(H) at H$sf ~ 2.73 kOe. In this study, M was measured using a vibrating 
sample magnetometer (VSM), in which the sample can be rotated so as that the in-plane 
magnetic field is applied at any skew angle tfj with respect to the easy axis. 
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B. Spin configuration for H parallel to the easy axis (ip = 0) 



For H applied in-plane parallel to the easy axis (if) = 0), the calculation of the ground- 



state magnetization profiles shows that, in a limited field range near H S sf = \j H e Ha + H\, 
the system admits two stationary configurations: the AF state, with collinear and antipar- 
allel layer magnetizations, and the surface spin-flop (SSF) state, with a non-uniform mag- 
netization profile characterized by a Bloch wall that nucleates near one of the film surfaces. 
In Fig. 1 we plot, in the neighborhood of Hssf, the field dependence of the reduced mag- 
netization, m(H) = M/M s = (1/N) J2iLi cos (f>i(H), and of the reduced energy e(H), given 
by Eq. ([I]). Using for the calculations H E = 9.80 kOe, H A = 0.98 kOe, N = 20 and if) = 0, 
the field of thermodynamic equivalence between the AF and the SSF state is found to be 
H th = H(eAF = cssf) = 3.02 kOe, while the boundaries of the metastability region are given 
by the fields H ini = 2.93 kOe and H sup = H S sf = 3.26 kOe. The AF state is the ground 
state for H < H th and is metastable for H th < H < H snp , while the SSF state is the ground 
state for H > iJ t h and is metastable for < H < H t h- The calculated value of the bulk 



SF field is H BSF = ^2H E H A + H\ = 4.50 kOe. 

In Figs. 2,a-d the calculated ground-state magnetization profiles are compared with those 
obtained from polarized neutron reflectivity measurements, as published in Ref. 
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at dif- 



ferent field values, ranging between and 5.5 kOe. The fields h = H/Hbsf were scaled with 
respect to the bulk SF field, where we set Hbsf — 4.50 kOe for the theoretical results and 
Hbsf ~ 4.14 kOe for the experimental ones. The orientation 4>AF(i) of the antiferromagnetic 
axis (i.e., the axis along which the magnetizations of two Fe adjacent layers are antiparallel) 
is plotted as a function of the Fe layer number for selected values of the applied field. Figure 
2e illustrates that <f>Ap(i) = \[4>(i) + <f>(i — 1)] — 90°, where <fi(i) is the angle formed by the 
magnetization Mj of the 2-th Fe layer with the field direction. As the measured spectrum 
extends only as far as the half order AF Bragg peak, which is determined by the antipar- 
allel components of the magnetizations, the orientation 0af(^) of the AF axis is obtained 
with more accuracy from the experiments than that of the individual layer orientations, for 
which the estimated error is up to 20°.— This representation clearly depicts the position and 
extent of the domain wall for the different fields and shows that agreement between theory 
and experiment is fairly good. 

From Figs. 2,a-d one sees that the basic features of the SSF transition in this film, char- 
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acterized by a low value of r « 1/10, are the following: i) for H > Hssf the deviations from 
the uniform AF spin configuration originate just at the surface layer whose magnetization is 
antiparallel to the field; ii) with increasing H, the surface-nucleated domain wall is pushed 
gradually into the middle of the film; Hi) for H > Hbsf a symmetric spin configuration is 
achieved, similar to the bulk SF one in the middle planes (while the spins at the surfaces, 
owing to the cuts of the exchange bonds, are less deviated from the field direction with re- 
spect to the bulk ones). Although not directly obvious from Fig. 2, the disco mmensuration 12 
at the center of the surface-nucleated domain wall, effectively dividing the AF order into two 
antiphase domains, which was strongly evidenced by the neutron data,— is also reproduced 
by the calculations. 

Note that no abrupt variations of the magnetization, except the one at Hssf, are found 
as the field intensity is increased, see Fig. 1: i.e., additional first-order transitions between 
C phases, intermediate between the AF and the SSF one,— ^ are not allowed by the low 
value of r = Ha/H e in this system. 

C. Spin configuration for H applied along an arbitrary direction (ip ^ 0) 

For a field applied in-plane along an arbitrary direction forming a skew angle if) ^ with 
the easy axis (i.e., both Hx and Hz are non zero), one has that M(H) = for H = 0. For 
< if) < 5°, M(H) increases slowly with increasing H, since for if) ^ the magnetizations 
of the two sublattices are no longer compensated. Upon further increasing H, a finite 
jump, signaling the onset of a first-order phase transition, was observed for a -^-dependent 
field value. The measurements were performed both upon increasing and decreasing H and 
showed a marked hysteresis for such small if) values: see Fig. 3. The magnetic susceptibility 
x(H) = dM/dH , obtained by numerically deriving the measured magnetization with respect 
to H, showed sharp peaks corresponding to the jumps in M(H): see Fig. 4, top. The 
full width at half maximum (FWHM, see later Fig. 5, top) of the measured susceptibility 
peak was found to be constant and very small (essentially determined by the instrumental 
resolution), signaling that for such small angle values (0 < if> < 5°) the SSF transition is of 
first order. 

As if} was gradually increased above 5°, M(H) became smoother (see Fig. 3), and the 
peak in x(H) (see Fig. 4, bottom) decreased in intensity, while the FWHM dramatically 
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increased on passing from 5° to 23°, as shown in Fig. 5, top. The latter feature strongly 
suggests that a crossover of the surface phase transition from first-order to second-order 
might take place for ip > 5°. 

In the light of this interpretation, it is however necessary to justify the persistence - up 
to the highest investigated value of ip (23°) - of a small hysteresis loop: see Fig. 5, bottom, 
where the measured dependence of the peak position of the magnetic susceptibility x{H) 
is shown as a function of the skew angle vp, both for increasing (full circles) and decreasing 
(open circles) magnetic field. 

To this aim, we observe that the coexistence of first- and second-order transition features 
was recently observed^ in single-crystal Lao.73Cao.27Mn03 perovskites exhibiting colossal 
magnetoresistance. The magnetization isotherms displayed a metamagnetic structure linked 
with a first-order transition, while field and temperature dependent ac susceptibility data 
presented a crossover line characteristic of a continuous transition.— 

In our case of a finite AF Fe/Cr(211) film, a similar effect, i.e. the coexistence of first- 
and second- order transition features, might be attributed to a distribution of values of the 
interlayer exchange, as well as of the anisotropy of the different layers in the stack, due to the 
presence of thickness fluctuations. In determining the observed small hysteresis loop, one 
cannot either rule out the role of defects (pinning centers) which inhibit the lateral motion 
of domains during the magnetization reversal process. 

In the following we will test if the experimental data can be explained in terms of a 
crossover from first- to second-order critical behavior by performing a theoretical calculation 
of the magnetization profile {<pi \ i — 1, • • ■ , N} for different values of H and ip ^ 0, using 
either of the two methods described in Sec. II. In the present case where r < 1, the map 
portrait is not chaotic and the two different methods give similar results. From the calculated 
spin configuration, one obtains m(H) = M/M s = (1/iV) Y^Li cos^j, x(H) = dM/dH and 
e(H), given by Eq. <Q. 

We find that for < ip < 4.75° the system admits two stationary configurations: a nearly 
AF state and a SSF state. The first one is stable only for H < H snp while the second one 
is stable only for H > H in {. By H t ^ we denote the field of thermodynamic equivalence at 
which the two states take the same energy. As ip increases, the width of the metastability 
region gradually reduces until, for ip > 4.75°, only one equilibrium configuration is found. 

In Fig. 6 the calculated susceptibility x{H) is shown for different values of the skew angle 
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if>. For if> < 3°, the peak in x{H) is a Dirac delta function, so the peak position is indicated 
by a vertical line. For ip > 5° the peak has a finite width and a finite height. As if> increases, 
the peak broadens and its height decreases. For clarity's sake, the peak position reported 
in Fig. 6 for if) < 3° corresponds to the calculated field of thermodynamic equivalence 
between the energies of the AF-like configuration and the non-uniform SSF configuration. 
It is just the position of this peak that is reported vs if) in Fig. 7 as the full-circle diagram; 
the other two diagrams plotted for ift < 4.75° represent the calculated field values H [n f (open 
squares) and H sup (open triangles) vs if). 

In the phase diagram of Fig. 7 one clearly observes that, upon increasing if>, the width 
of the metastability region gradually shrinks until, for if> greater than a critical value 
V'max (film) =4 . 75° , the film admits only one equilibrium state. Thus, for if) > -?/> max (film) we 
expect the SSF transition to become continuous. The calculated value of ^ max (film)=4.75° 
turns out to be in remarkable agreement with the value ~ 5° estimated from the experi- 
mental results on the basis of the strong increase observed in the FWHM of x{H). This 
fact provides support for the hypothesis of a crossover from first- to second-order critical 
behavior for the SSF transition in a skew field. 

The calculated value of ^ max (film)=4.75° should also be compared with its bulk counter- 
part. In the bulk case, the field-induced phase transition of a uniaxial antiferromagnet in the 
presence of a skew field forming an angle if) with the easy axis was theoretically studied by 
Rohrer and Thomas. 30 They predicted the first-order bulk SF transition to become continu- 
ous for if) > ■^ max (bulk)=tan~ 1 [_f/' J 4/(2if£; — Ha)]. However, in MnF 2 , where r m 1/100, the 
critical angle turns out to be as small as ~ 0.4°. In the case of the Fe/Cr superlattice under 
study, the ratio r is nearly an order of magnitude higher than in MnF 2 , so that an appre- 
ciable critical angle ?/> max (bulk)^ 3° is estimated. The calculated value of ^ max (film) =4.75° 
is nearly twice the bulk value. This can be qualitatively understood considering that, for 
not too high values of r, in the bulk the critical angle is essentially determined by the ratio 
Ha/ (2He), while in the film the effective exchange field at the surface is halved with respect 
to the bulk. 
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IV. DISCUSSION 



In this work we have investigated the transition, induced by a magnetic field H with 
arbitrary direction, between the antiferromagnetic phase and the surface spin-flop phase of 
an epitaxial Fe/Cr(211) superlattice with t Fe = 14 A, t Cr = 11 A and N = 20 repetitions. 
The system is characterized by a rather small value (r pa 1/10) of the ratio r = Ha/ He 
between the uniaxial anisotropy field Ha and the exchange field He, yet much greater than 
the value (r pa 1/100) pertinent to usual bulk antiferromagnets like MnF2 and C^Os. 

For an external field applied parallel to the easy in-plane axis (ip = 0), the layer-by- 
layer spin configurations measured by polarized neutron reflectometry were found to be in 
remarkable agreement with theoretical calculations, performed in the framework of a mean- 
field ID model of the superlattice stack. 

For a field applied in-plane along an arbitrary direction forming a skew angle if) ^ with 
the easy axis, the superlattice magnetization M(H) was measured using magnetometry. The 
phase diagram of the film was calculated in order to check the possibility of a crossover of the 
surface phase transition from first- to second-order to take place for -0 max (film), similarly to 
what was predicted decades ago by Rohrer and Thomas^ for an AF bulk system in a skew 
field. Indeed we calculated ^ max (film)^ 4.75°, to be compared with the experimental value, 
ip pa 5°, at which the jump in M(H) starts smoothening and the FWHM of the measured 
magnetic susceptibility displays a dramatic increase with increasing if>. Owing to the cut 
of exchange bonds at the film surfaces, the calculated value of ip maiX (fi\m) turns out to be 
nearly twice its bulk counterpart, ?/' max (bulk)^ 3° M The latter value is much higher than 
the ones predicted for ordinary bulk antiferromagnets {e.g., ^ max (bulk)ps 0.4° for MnF 2 and 
pa 0.015° for Cr 2 3 ). 2 ° 

From the comparison between our experimental and theoretical results we conclude that 
a crossover between first- and second-order critical behavior is easier to be observed by 
magnetization measurements in Fe/Cr superlattices, thanks to the much higher value of the 
ratio r = Ha/ He between the anisotropy and the exchange fields in such an artificially 
grown system with respect to ordinary bulk antiferromagnets. 

The interpretation of the experimental data proposed above needs, however, further in- 
vestigation to be conclusive. In particular, a quantification of the magnetic domain structure 
in the presence of a structurally rough interface is required. While our previously published^ 
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polarized neutron reflectometry data proved unambiguously that only one type of domains 
is present in Fe/Cr(211) superlattices with an interfacial roughness of ~ 4 A for zero field, 
such a clear-cut evidence is lacking in the case of a nonzero field, applied in plane along an 
arbitrary direction. In conclusion, while we do not claim a quantitative accuracy for our 
theoretical results, nevertheless we believe that the main features of the spin-flop transition 
in the Fe/Cr superlattice have been captured by our "ideal" model (i.e., characterized by 
structurally smooth and uniformly magnetized layers). 
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FIG. 1: (Color online) Calculated field dependence of the reduced magnetization, m(H) = M/M S: 
and reduced energy e(H), given by Eq. (1), of an antiferromagnetic film with N = 20 planes, 
for H in the neighborhood of the surface spin-flop transition. The field is applied along the easy 
axis (■0 = 0); the exchange field and the anisotropy field are, respectively, He = 9.80 kOe and 
Ha = 0.98 kOe. The thick (thin) line refers to increasing (decreasing) magnetic field. The field 
of thermodynamic equivalence between the energies of the collinear AF and the non-uniform SSF 
configuration is H t h = 3.02 kOe. The collinear AF state is stable for H < H sup = 3.26 kOe; the 
non-uniform SSF state is stable for H > H m { = 2.93 kOe; they have the same energy at H t h- 
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FIG. 2: (Color online) a-d) Comparison between the calculated layer-by layer configuration (solid 
lines) and the experimental ones (markers and dashed lines) deduced from polarized neutron reflec- 
tometry data in a [Fe(14 A)/Cr(ll A] x 2o superlattice. The orientation (f>AF of the antiferromagnetic 
(AF) axis (i.e., the axis along which the magnetizations of two adjacent Fe layers are antiparallel) 
is plotted as a function of the Fe layer number for selected values of the reduced applied field 
(h = H/Hbsf); e) Schematic drawing showing that 0af(*) = 5 + <H^ ~~ 1)1 ~~ 90°, where <j>(i) 
is the angle formed by the magnetization Mj of the i-th Fe layer with the field direction. The 
magnetic field is applied parallel to the easy axis (ip = 0). 
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FIG. 3: Experimental VSM data for the magnetization M{H) of a [Fe(14 A)/Cr(ll A)] x20 super- 
lattice vs H, applied along a direction which forms a skew angle Y> with the easy in-plane axis. 
The two different curves refer to ip = 1° and ip = 23°, respectively. 
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FIG. 4: (Color online) Experimental data for the magnetic susceptibility x{H) = dM/dH of a 
[Fe(14 A)/Cr(ll A)] x 2o superlattice, obtained by numerically deriving the VSM data in Fig. 3, 
both for increasing field (thick line) and decreasing field (thin line). Top diagram: = 1° ; bottom 
diagram: tp = 23°. In each diagram, the arrow at the lower (higher) field denotes the peak position 
for decreasing (increasing) H. 
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FIG. 5: Top diagram: experimental peak intensity (open squares) and FWHM (full squares) of 
the magnetic susceptibility x(-fT) °f a [Fe(14 A)/Cr(ll A)] x 2o superlattice vs the skew angle ip 
formed by the applied magnetic field with the easy in-plane axis. The lines are guides to the eye. 
A dramatic increase in the FWHM is observed for > 5°. Bottom diagram: experimental peak 
position of x{H) vs if), both for increasing (full circles) and decreasing (open circles) field intensity. 
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FIG. 6: (Color online) Calculated X (H) = dM/dH of an AF film with iV = 20 and H E = 9.80 
kOe, Ha = 0.98 kOe vs H applied in-plane along an arbitrary direction. The different curves refer 
to different values of the skew angle V formed by the external magnetic field with the easy axis. 
For ip < 3°, the position of the reported peak of x{H) is indicated by a vertical dashed line and 
corresponds to the field of thermodynamic equivalence i? t h between the energies of the AF-like and 
the non-uniform SSF configuration. 
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FIG. 7: (Color online) Calculated phase diagram of an AF film with N = 20 and He = 9.80 kOe, 
Ha = 0.98 kOe. The peak positions of x{H) = dM/dH, scaled with respect to Hssf, are reported 
vs the skew angle tp formed by the external magnetic field with the easy axis. For tp > 4.75° 
the SSF transition is predicted to become continuous. The calculated field of thermodynamic 
equivalence H t ^ is reported vs if) as the full-circle diagram; the other two diagrams plotted for 
i\) < 4.75° represent the calculated field values H m f (open squares) and H sup (open triangles). 
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